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Abstract 

A potential model to describe the total cross section of nuclear fusion reactions at low energies is 
proposed. It is shown that within an approach with a simple, single barrier potential a satisfactory 
description of data is hindered, while with models with more complex interactions, e.g., with two- 
barrier potentials, a good description of data can be achieved in a large interval of colliding nuclei 
and center of mass energies. In particular, the two-barrier model allows to describe data at low, 
sub-barrier energies where data exhibits a steep falloff of the cross section. It is also shown that 
the position of the two barriers is almost independent up on masses of colliding nuclei. Comparison 
with available experimental data is presented as well. 
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I. INTRODUCTION 



Investigation of fusion reactions at low center of mass energies (i? ^ 30 — lOOMeV) is 
of interests not only in context of study of the mechanism of nucleus-nucleus interaction, 
but also supplies unique information on astrophysical processes on stars, in particular on 
the so-called astrophysical S-factor, which is directly related with the e.g., flux of solar 
neutrinos. Experimental measurements at energies well above the Coulumb barrier started 
a few decades ago [l, 2] and shown that the cross section exhibits a smooth behaviour as a 
function of the center of mass energy. However, measurements of very small cross sections 
{aER ~ lO^^mfe), near and below the Coulomb barrier js-Q] have shown an unexpected steep 
falloff of the cross section, named as "hindrance" effect. Theoretical efforts to describe such 
a shape generated a number of models 9l-|l5l| within which different reaction mechanisms 
have been suggested to reconcile data at extreme sub-barrier energies. Since different models 
describe data equally well, in order to distinguish the proposed mechanisms one needs more 
measurements at even lower energies; a comparison with model predictions will allow to 
enlighten the peculiarities of the reaction mechanism. 

In the present paper a further development of the so-called " critical distance model" 

m 

is presented. Remaind, that within such a model, in order to create the compound nucleus, 
the colliding nuclei have to approach some "critical distance" Rcr with = ^cri^Y^ + ^2^^) 
{r^r ~ 1/^)- At such distances the colliding nuclei loss their individual structure forming 
a single highly excited system. Since at the initial stage the collision is mainly superficial 
the deformation nuclei can be neglected and, correspondingly, the ion-ion potential can be 
calculated within the so-called "frozen approximation" S]- It should be noted that when 
the density in the overlap region becomes larger than equilibrium density, the Pauli principle 
forbids further penetration of nuclei in to each other. This gives rise of a repulsion part in the 
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ion-ion potential [16|, 117| at such critical distances. At smaller distances r < R^r the colliding 
nuclei forming the compound system become deformed so that the densities cannot be 
considered "frozen" anymore. This leads to a drastic attenuation of the potential up to values 
Qr = {McN — Ml — M2) c^, where Qr is the threshold energy of the reaction and M denotes 
the corresponding mass of the nuclear system. It is clear, that if the critical distance Rcr is 
smaller than typical positions of the Coulomb barrier, this results in a two barrier potential; 
besides the Coulomb barrier here appears another, inner barrier Vc{Rcr)- The Coulomb 
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barrier dominates at small values of atomic numbers however with increasing of Z the 
first, inner barrier becomes similar and even larger than the Coulomb barrier. Usually the 
strength of potentials Vc{Rcr) is characterized by a parameter ?7cr. {j] = Zi Z2/ (aY'^aI^^)) [l^ 
which demonstrates that already at r/cr > 71.3 the first barrier, Vc{Rcr), becomes larger than 
the Coulomb one, Vb{Rb)- For light and intermediate nuclei Vc{Rcr) is always smaller than 
Vb{Rb) and it contributes only at large angular momenta and high energies. 

The aim of this work is to investigate the effects of the two barriers in the ion-ion potential 
and to explain the observed hindrance as effects of resonance tunneling through such a two- 



humped potential 



lOj. It will be shown that the inner barrier manifests also in deep sub- 



barrier cross sections of medium magic and semi- magic nuclei even at 77 < rjcr- Particular 
attention is paid to fusion reactions of nuclei with mass number closed to the magic ones. 
In this case according to the nuclear shell model the single particle energy levels are closed 
and the nucleus acts as a hard sphere impeding the deformation during the fusion process. 
As a result, the tunneling of undeformed ions occurs without forming a compound system 
and the interference between two barriers can become important. This effect is studied in 
the present paper in details. 

Since the procedure of finding the exact solution of the Schrodinger equation (actually, 
this is a system of coupled equations) is a cumbersome and rather tedious task, in the present 
approach a semiclassical approximation is adopted within which the system is considered 
uncoupled and each equation solved separately. The resulting solution contains a number of 
free parameters which determine the general form of the inner potential. These parameters 
are to be found from a combined fit of the experimental data for different pairs of colliding 
nuclei. The possibility to apply the proposed model to other kind of processes with heavy 
ions is discussed as well. 



II. FORMALISM 



The equations for radial wave functions Raiai''^)-, describing the relative motion of the 



two ions in the coupled channels method can be written as 



19| 



RcLAr) = '^j:Vap{r)RpL,{r). (1) 
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Here La is the orbital angular moment, V^(r) is the diagonal part of ion- ion potential 
in the a channel with excitation energy Ea and Vapir) is the coupling potential between 
channels. Eventually, /i is the reduced mass of the system. 

As already mentioned, the coupled system ([T]) can be decoupled and solved by considering 
the semiclassical approximation 20|]. Within this approximation the system is spinless for 
both channels and the energy levels are considered as a two-level degenerated system in the 
ground state. Introducing the notation x^l\^) = Roiir) ± -RiL(r) and using the equalities 
Vo('") = Vi{r) and Voi(r) = Vio{r), one obtains 



+ -fc — 



{Vo{r)±Voi{r)) 



(2) 



2/i dr"^ 2jj,r'^ 

In the simplest case when the fusion process is supposed to be defied as a conserved flux 
of ions tunneling throughout the Coulomb barrier the system can be solved analytically 
resulting in the known Wong formulae 2l|]. Then in this case the fusion cross section can 
be expressed as 
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afUE) = {ln[l + exp '^{E - V,{Rb) - V,,{Rb))]] + 



AE 
huR^ 



B 



2ti 



^ln[l + exp— (E-yo(i?B) + K)i(i?B))] 



(3) 



In on]) the radii in both channels have been taken the same. 

In solving the equation ([2]) the quasi-classical method WentzelKramersBrillouin (WKB), 
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23|, has been used. Accordingly, the 



originally proposed to describe fission reaction 
penetrability T can be written 

T = TcTb/{{1 + [(1 - Tc){l - Tb)]'/^Y cos' u + {1 - [{1 - Tc)il - Tb)]'/'}' sin^ u), (4) 

where for the sake of brevity in the above equation the angular momentum L dependence 
has been suppressed; the subscripts B and C stand for the Coulomb and inner barriers, 
respectively. The corresponding penetrability Tc,b within the WKB-approximation at sub- 
barrier energies is 

Tc,B {E) = {1 + exp [2 ScAE) , (5) 

where the actions Sc,b{E) and phases u are defined by integrals over the interval between 
corresponding return points in ion-ion potentials VL{r) 



ScAE) = J\2fi[VL{r)-E]/h' 



1/2 



dr 



(6) 



dr. 



(7) 



p = j \2ii[E-VL{r)]/h^\"^ 

At larger penetrability, Tc^b > 0.45, the tunneling probabilities can be parametrized 
analytically by known Hill- Wheeler formulae with near-barrier parameters of curvatures 



C,B 



dr"^ 



1/2 

r=Rc,B 



(8) 



TcAE) = 1 + exp 
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-(E-Vcb) 



(9) 



Here Vc\b, Rc,b are potentials and radii of inner and Coulomb barriers of each partial wave. 
Then the averaged, over z/ in eq. (jlj), penetrability can be written as 



T = TcTB/[l-{l-Tc)il-TB)]. 



(10) 



Exactly the same expression for the penetrability can be obtained by summing ingoing and 
outgoing fluxes in the interval between the two barriers. Note, that by averaging in obtaining 
ffTOj) all the interference effects have been lost. So that a comparison of results by unaveraged 
formula (jl]) and by equation (fTOl) will reflect the role of interference effects in reactions of 
heavy ion fusion at low energies. 



III. DATA ANALYSIS 

We have analyzed fusion reaction for a variety of combinations of colliding ions available 
for experimental studies. The corresponding reactions together with some their specific 
characteristics (Q-reaction, minimal energy, Bass potential and angular momentum of the 
evaporation residue) are listed in Table [H In our approach these quantities serve as input 
parameters. Our analysis consists on changing the form of the of ion-ion potential and 
computing the fusion cross section with further fit of the potential parameters to obtain a 
best description of data. As the outer part of the Coulomb potential the Wood-Saxon form 
with parameters V, R, a has been used. The coupled channel potential has been taken in the 
form 



dVoir) 
dr 



+ 



Z,Z,e'R',-'/r'^' 



2A + 1 



(11) 
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TABLE I: The considered reactions and the relevant parameters used in the numerical calculations. 
Qr is the threshold energy of the reaction, Emin is the minimal energy at which the experimental 
measurements have been performed, Jer is the parameter used to distinguish the two deexcitation 
process, eq. (fTSl) . Vsass is the Bass potential used in most phenomenological calculations. 



Reactions 


-QniMeV) 


VsassiMeV) 


EminiMeV) 


Jer 


Refs. 


64 +64^^ 


48.71 


95.7 


85.55 




m 


60^^ +89 y 


90.64 


129.9 


121.4 








154.0 


180.9 


169.6 


17.0 


[4] 


28 Si +64 


1.78 


52.3 


44.0 




M 




46.49 


76.0 


62.7 


6.0 





where A is the multipolarity of the level and the dynamical deformation parameter /3q, is 
taken from independent experiments. In the present calculations the quantity {(3aRo) / V^ti 
has been taken constant and the same for all reactions listed in Table [T] {fiaRo/ = 0.3). 

The form of the inner part of the central potential imitating the repulsive interaction is 
determined, in the present paper, by the expression 



AVo{r) = ViiR,^-r)^. (12) 

which, at r < R^ {Rm is determined by the position of the minimum of the central potential) 

16| . To compensate the increase of the potential (fT2|) at low 



is similar to one used in Ref. 



distances (increase of R^ — t), an additional damping factor 



^(r) = ( 1 + exp ^— ^V' (13) 



A 

has been introduced with parameters Ri A as free ones. In such a way the total potential 
results in a continue function of the distance r with continue first derivatives, as it should 
be for a nuclear potential. Explicitly, the potential reads as 



Voir) ± K)i(r) = Qr ■ g{r) + (K)(r) ± K)i(r)) ■ (1 - g{r)). (14) 

The potential (fT4|) is used to calculate the nuclear penetrability T, eq. Besides 
penetrability, the fusion cross section is determined by the probability of deexcitation of the 



TABLE II: The final values of the phenomenological parameters of the model obtained by fitting 
the experimental data. The quantities V, R and tq are the magnitude of the barrier, its position 
and nuclear parameter for Coulomb and inner potentials, respectively. 



Reactions 


Vb 


Rb 


^0 


hUB 


Vc 


Rc 




hujc 


64 +64 


93.55 


11.28 


1.41 


3.57 


85.11 


6.68 


0.84 


12.17 


60^^ +89 y 


129.65 


11.50 


1.37 


4.33 


122.33 


6.92 


0.83 


9.76 




173.63 


12.35 


1.37 


3.46 


178.23 


7.87 


0.87 


10.56 


28 5i +64^^ 


51.26 


10.35 


1.47 


3.92 


45.05 


6.04 


0.86 


19.7 


160+208P6 


74.47 


11.98 


1.42 


4.76 


61.86 


7.44 


0.88 


14.83 



compound system via fission and evaporation processes. In a compound system formed by 
light nuclei the evaporation process prevails, while for heavy components the fission reaction 
dominates. In the medium region both processes are of an equal importance. To distinguish 
between these two processes of deexcitation we use a simplest model, suggested in j25| . 
It consists in merely multiplication of the penetration probabilities in(jl]) and ([9]) by the 
"deexcitation" factor 

Wer{L) = (l + exp ^^^=^) ' , (15) 

where angular momentum Jer separates fission and evaporation residue cross sections, cf. 
Table E 

Having determined the form of the cross section we fitted the free parameters to obtain 
a good description of data in a large region of energies and nuclear masses. The results of 
our analysis are shown in Table [TTl where the fitted values of the relevant parameters of the 
model are listed. The main result from the Table can be formulated as follow: since for all 
combinations of the colliding ions the difference Rb — Rc is almost constant (~ 4.5 fm) the 
distnace between two barriers is approximatively the same for all the considered reactions. 
Moreover, as expected Rb — Rc is roughly twice the diffuse edge of nuclei. 

As an example, in Figs. [T] we present results of the cross section calculated (solid lines) 
within our model with parameters listed in Table [TTl for two combinations of colliding nuclei. 
Since in such kind of reactions the colliding nuclei are considered as light, the theoretical cross 
section is compared with the experimentally measured evaporation residue cross sections. It 
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FIG. 1: The calculated cross sections Eq. ^ for the reactions ^^Ni Ni, left panel and for the 
reaction ^^Ni Y, right panel. Solid curves correspond to full calculations with including the 
interference effects in penetrability, eq. (SI); dashed curves are the results of calculations without 

un 

interferences, eq. (jlOp . Experimental data correspond to the measured evaporation residues [3, l5|]. 

can be seen that within two-barrier model a good description of data can be achieved in the 
whole interval of considered energies for both combinations of colliding ions. It should be 
pointed out that the hindrance effect has been obtained only by an accurate account of the 
interference effects, as above discussed. Without interference, i.e. cross section by formula 
(fTOj) (dashed lines in Figs. [1]) an increase of the slope of the cross section at low energies 
cannot be obtained. Our results are in an agreement with the ones of Ref . |t2^] , the fusion 
cross sections for the reactions ^^Zr +^'^^Sn, ^^Kr Sh +'^^^Pb, where rj > rjcr have 
been analyzed. 

As an illustration of the obtained potential in Fig. |2]we plot the fusion potential without 
the Coulomb barrier, V{r) — Vb as a function of the distance from the Coulomb barrier 
r — Rb for five combinations of the colliding particles. It is clearly seen that the position 
of the barriers relative to Rb is roughly the same. This is an important result since if so, 
predictions for other sets of colliding nuclei can be easily performed. 

Another important result obtained during the phenomenological fit is that at = Qb. 
the phase v in eq. ([7]) becomes ~ 7r/2, which is a clear evidence of the fact that in the 
inverse to the fusion channel the decaying system was in its first quasi-stationary state. It 
means that the obtained potential ((Hj) can be successfully applied in describing the decay 
of heavy nuclei into heavy clusters. 
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FIG. 2: The phenomenological potential, eq. (jl2p . as a function of the "scahng" variable r — Rb 
for a set of colliding ions. The scaling- like behaviour of the inner potential is clearly seen at 
r — Rb — —4:.5fm. 

IV. CONCLUSION 

The fusion processes of two nuclei have been analyzed within a model with complex 
shaped ion-ion potential. It has been shown that the observed hindrance effect in the fusion 
cross section at extremely low energies can be described by a two-barrier potential. The 
corresponding shape and values of the phenomenological parameters have been obtained by 
a fit of experimental data and it has been shown that the potential depends rather on the 
relative distance from the Coulomb barrier that the distance between the colliding ions. It 
has been found that the interference contribution in to the penetrability T plays a crucial 
role in understanding the drastic change of the slope of the cross section at extremely sub- 
barrier energies. The obtained phenomenological potential possesses a scaling-like behaviour, 
in sense that it is almost the same for different ions at the same values of the distance from 
the Coulomb barrier. This allows to apply our model to other kind of reactions, e.g. the 
heavy cluster decay processes. 
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